Accepted for Publication in The Astrophysical Journal 
Preprint typeset using LATf^X style emulateapj v. 11/12/01 



SPITZER IRAC SECONDARY ECLIPSE PHOTOMETRY OF THE TRANSITING EXTRASOLAR 

PLANET HAT-P-1B 



O 
O 

(N 

> 
O 



CM 

6 



(N 
> 
oo 

(N 
(N 



o 



13 



Kamen Todorov 
William C 



1,2,3 

Bowman 4 , 



, Drake Deming 2 , Jospeph Harrington 4 , Kevin B. Stevenson 4 , 
Sarah Nymeyer 4 , Jonathan J. Fortney' 1 , and Caspar A. Bakos 6 ' 7 

Accepted for Publication in The Astrophysical Journal 

ABSTRACT 

We report Spitzer/IRAC photometry of the transiting giant exoplanet HAT-P-lb during its secondary 
eclipse. This planet lies near the postulated boundary between the pM and pL-class of hot Jupiters, 
and is important as a test of models for temperature inversions in hot Jupiter atmospheres. We de- 
rive eclipse depths for HAT-P-lb, in units of the stellar flux, that are: 0.080% ± 0.008% [3.6/zm], 
0.135% ± 0.022% [4.5/Ltm], 0.203% ± 0.031% [5.8/im], and 0.238% ± 0.040% [8.0/mi]. These values are 
best fit using an atmosphere with a modest temperature inversion, intermediate between the archetype 
inverted atmosphere (HD 209458b) and a model without an inversion. The observations also suggest 
that this planet is radiating a large fraction of the available stellar irradiance on its dayside, with little 
available for redistribution by circulation. This planet has sometimes been speculated to be inflated 
by tidal dissipation, based on its large radius in discovery observations, and on a non-zero orbital ec- 
centricity allowed by the radial velocity data. The timing of the secondary eclipse is very sensitive to 
orbital eccentricity, and we find that the central phase of the eclipse is 0.4999 ± 0.0005. The difference 
between the expected and observed phase indicates that the orbit is close to circular, with a 3<r limit of 
e cos uo\ < 0.002. 

Subject headings: stars: planetary systems - eclipses - techniques: photometric 



1. INTRODUCTION 

Infrared light from transiting extrasolar planets can be 
measured using high precision photometry when the plan- 
ets pass behind their stars. Most of these secondary 
eclipse me asurements to date have u s ed the Spitzer Space 
Telescope (ICharbonneau et al.l 120051: iDeming et all [2005, 
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Knutson et al.ll2008ll2009h . Spitzer data in th e four bands 



of the IRAC instrument (|Fazio et alj I2004D define the 
planet's spectral energy distribution over the 3.6-8.0 /im 
wavelength region, where water vapor plays a principal 
role in shaping the spectrum. Spitzer studies have shown 
that there are at least two classes of close-in extraso- 
lar giant planets, differentiated by th e temperature gradi- 
ent with height in their atmosphere (IHubenv et al.l 120031 : 
iFortnev et all2006l . 120081 IBurrows et al.ll2007l . l200l . One 
class (pM) exhibits a tempe rature inversion at h i gh al- 
titude in their atmo sphere (jKnutson et al.l 120081 , [2009; 
IMachalek et~aLl 12008th The inversion affects the emer- 
gent spectral energy distribution by causing the water 
bands to appear in emission. The cause of the inver- 
sion is believed to be absorption of strong stellar irradi- 
ation in the visiblej poss i bly by TiO /VO (jHubenv et alj 
20031 IFortnev et aT]|2006L [20081) . but other absorbers re- 



m ain possible (IBurrows et al J 12001 IZahnle et~alll2009f) . 

IMachalek et all (|2008t ) found that XO-lb exhibits an 
inverted atmospheric structure, in spite of being nomi- 
nally below the stellar irradiance level projected to de- 
fine the transition between the pM and the (non-inverted) 
pL classes. The pM/ pL transition corres ponds to the 
condensation o f TiO (jFortnev et alj 120081 ). Recently, 
IMachalek et al.l (|2009h found that XO-2b exhibits a weak 
temperature inversion, consistent with being near the 

M/pL transition . Like XO-lb and XO-2b, HAT-P-lb 
Bakos et al.ll2007[ ) lies at the lower edge of that predicted 
boundary, so it provides an important additional test of 
the propensity toward inversion at that level of irradia- 
tion. HAT- P-lb may also exhib it a non-zero eccentricity 
to its orbit ([Johnson et al. 2008), that may cause internal 
energy generation via the dissipation of tidal stress. 

In this paper, we report Spitzer photometry of HAT- 
P-lb's secondary eclipse in the four IRAC bands. We 
use these data to investigate the atmospheric temperature 
structure and energy balance of this planet, and to place 
more stringent limits on the ecc entricity of its orbit via th e 
timing of the secondary eclipse ([Charbonneau et al. 2005). 

2. OBSERVATIONS 

We observed HAT-P-lb during two secondary eclipses 
using Spitzer/IRAC. On 2006 December 28, we observed 
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an eclipse at 4.5 and 8.0 /jm for 356 minutes. To avoid 
saturation for this relatively bright star, each frame com- 
prised two 2-second exposures in stellar mode at 4.5 /im 
(3212 exposures total), simultaneously with one 12-second 
exposure at 8 /im (1606 exposures total). On 2007 Decem- 
ber 29, we observed an eclipse at 3.6 and 5.8 /im for 333 
minutes (1510 total exposures at each wavelength). For 
that eclipse, we elected to increase the efficiency at 3.6 /im 
and avoid saturation by placing the star at the corner of a 
pixel, thus spreading the light over 4 pixels. We used 12- 
second exposures at both wavelengths, in full array mode. 
For both eclipses, we centered the star on a portion of 
the array chosen to avoid known bad pixels and scattered 
light from bright stars imaged onto other regions of the 
focal plane. 

3. PHOTOMETRY 

All of our photometry used version S18.5.0 of the Basic 
Calibrated Data from the Spitzer pipeline. We applied the 
recommended factors to correct for the variation in flat- 
field response to point sources versus extended sources, 
and the variation in pixel solid angles, as described in Sees. 
5.3 and 5.6.2, respectively, of the IRAC Data Handbook 
V3.0. We ran our analysis both with and without these 
corrections, and we find that they change our final results 
by less than the la error in the eclipse depth at all wave- 
lengths. Their largest effect is at 4.5 /xm, where they in- 
crease the eclipse depth by 6 x 10~ 5 in units of the stellar 
flux. All of the results quoted in this paper include these 
corrections. 

To facilitate the error analysis, we converted the inten- 
sities in the images to electron numbers, using the calibra- 
tion information in the FITS headers. We corrected en- 
ergetic particle hits by comparing each pixel to a median- 
filtered time series of that pixel's intensity, using a 5-frame 
resolution. We replaced individual values exceeding the 
median by more than 4<r with the median value. Some en- 
ergetic particle hits that overlie the stellar image were not 
well corrected, and resulted in outlying intensities in the 
photometric time series; those values were omitted when 
fitting eclipse curves to the data. 

All of our results are based on aperture photometry, but 
the details differ with wavelength, as described below. At 
all wavelengths, we centered the aperture on each stellar 
image by fitting a parabola to the brightest three points 
in the stellar profile. The image was summed in X to de- 
fine the profile as a function of Y, and vice versa. We also 
found the center of each image by fitting a 2-dimensional 
Gaussian to the core of the image point spread function, 
but this did not improve the results over the parabolic fit. 
We varied the size of the photometry aperture, and used 
the value that minimized the scatter in the time series. 

At each wavelength, we subtracted the background due 
to solar system zodiacal thermal emission. Since HAT-P-1 
has a brighter companion star 11 arcsec distant, there was 
also a contribution due to diffracted light from the com- 
panion. Based on model PRFs provided by the Spitzer 
Science Center, we expect a diffracted contribution of 
~ 0.3MJy/sr at the position of HAT-P-1. This contri- 
bution is only weakly dependent on wavelength, because 
the decreasing brightness of the star at longer wavelength 
works in the opposite direction to the effect of diffraction. 



However, the zodiacal background is significantly wave- 
length dependent, so diffracted light from the companion 
star can dominate at 3.6 and 4.5 /im, decreasing to ~ 15% 
of the zodiacal background at 8 /im. We therefore mea- 
sured the background at a symmetric position on the op- 
posite side of the companion star from HAT-P-1, using an 
aperture of the same size and shape. Since this relatively 
small aperture (typically about 3 pixels in radius) encom- 
passes relatively few electrons from the background, we 
increase the precision of the background measurement by 
fitting a parabola to the background time series, and using 
the value of that fit for each frame. Apart from random 
error, the background variation was quite gradual and was 
well represented by the parabolic fit. 

3.1. 3.6 and 4-5/im 

Spitzer photometry at 3.6 and 4.5 /im is known to be 
affected by pixel position, wherein the value from aper- 
ture photometry is a functi on of the location of the stellar 
centroid within the pixe l (|Morales-Calderon et alj [2006; 
ICharbonneau et aT1l2005D . While this is true for our pho- 
tometry at these wavelengths, an additional factor is im- 
portant for our 3.6 /im observations. The upper panel of 
Figure 1 shows the photometric intensity of the star at 
3.6 /im, as a function of the Y-position (the X-position 
effects were not significant). The bulk of the data were 
collected with the star displaced significantly from pixel 
center, as was our intention when planning the observa- 
tions. The median value for the maximum per-pixel in- 
tensity in these images is about 130,000 electrons, approxi- 
mately equal to the 1% non-linearity point for this detector 
(see the Spitzer Observing Manual, Sec. 6.1.3). However, 
due to pointing jitter the star sometimes wanders closer to 
pixel center, and the intensity drops as saturation begins 
to affect the data (red points on Figure 1). Those points 
were not included in our analysis. The remaining data 
still exhibit variations in intensity that are dependent on 
pixel position, but the variation is complex because of the 
location of the star near pixel boundaries. We required a 
sixth-order polynominal in Y-position to fit this relation; 
incorporating the X-position in the fit did not improve it. 
We subtracted the fit (blue line on Figure 1) from the un- 
saturated photometric points, to decorrelate the variation 
due to pixel position. 

The lower panel of Figure 1 plots the uncorrected in- 
tensity versus orbital phase, showing that the eclipse is 
visible even without the correction for pixel position. This 
panel also shows that the omitted points (in red) occur 
periodically as the oscillation in telescope pointing moves 
the stellar centroid toward the center of the pixel. Since 
most of the correlation between intensity and pixel posi- 
tion depends on the Y-coordinate of the star, we used a 
square photometric aperture at both 3.6 and 4.5 /im, in- 
cluding fractional pixels at the edge. We reasoned that 
a square aperture would provide the cleanest isolation of 
pixel-position effects that are predominately dependent on 
the Y-coordinate. We found that the photometric scatter 
at 3.6 /*m was minimized using an aperture 5 pixels on a 
side (2.5 pixels in 'radius'). 

Our photometry at 4.5 /im shows a much weaker pixel 
position effect (not illustrated) than at 3.6 /im, but 
does not involve saturation. The 4.5 /im data show 
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a normal dependence on d istance from pixel center 
(jMorales-Calderon et al.|[2006h . but again with a stronger 
dependence on the Y- than the X-coordinate. The correla- 
tion coefficient between the image <5Y-position and inten- 
sity is -0.3, indicating a weak correlation, but statistically 
significant considering the 1604 degrees of freedom. We 
removed this correlation with a linear relation between in- 
tensity and radial distance, and the minimum scatter was 
achieved using a box 7 pixels on a side (3.5 pixels 'radius'). 
We explored other methods to achieve this decorrelation, 
such as using a function of Y only, and both X and Y 
separately, but the minimum \ 2 was achieved using radial 
distance. Figure 2 shows the 3.6 and 4.5 /im photometry 
plotted as a function of orbital phase, with the decorre- 
lation functions overplotted. Note that the decorrelation 
functions are smooth when plotted versus pixel position, 
but they become more jagged when plotted versus orbital 
phase as in Figure 2. 

Both the 3.6 and 4.5 photometry achieve a precision 
near the photon limit after the decorrelations. Specifi- 
cally, the scatter at 3.6 /im is 0.00131, which is merely 
6% greater than the photon-limited value. At 4.5 /im, the 
scatter of 0.00533 exceeds the photon-limit by 11%. The 
eclipse fitting procedure at every wavelength always used 
the unbinned data, but some plots show binned data for 
clarity. Figure 3 shows binned photometry for both 3.6 
and 4.5 /im, with the best-fit eclipse curves and ±1ct error 
limits. 

3.2. 5.8 and 8.0[im 

Photometry at 5.8 and 8.0 /im did not exhibit a de- 
tectable pixel position effect, there being no significant 
correlation between intensity and either the X- or Y- 
coordinate of the image. We performed the photometry 
using both square and circular apertures. Although our 
results do not depend significantly on the adopted shape 
of the photometric aperture, we elected to use a circular 
aperture at both 5.8 and 8.0 /im. We found a minimum 
scatter in the photometry using aperture radii of 2.4 and 
2.8 pixels at 5.8 and 8.0 /im, respectively. 

3.3. Eclipse Amplitudes 

Following the background subtraction and aperture pho- 
tometry, we divide the time series at each wavelength by 
its average value. This places the results in contrast units, 
i.e., relative to the stellar flux. 

We generated an eclipse curve nu merically, using the 
stell ar and planeta r y par ameters from lWinn et al.l (|2007| ) 
and I Johnson et aTl (|2008l ). The numerical code was tested 
for trans it curves by comp a ring t o the analytic formulae 
given by iMandel and Ago! (|2002t ), and it attains an ac- 
curacy of 10~ 6 , more than sufficient for our purpose (an 
eclipse curve has the same shape as a transit curve in the 
limit of zero limb darkening). In fitting to each set of 
photometry, we scale the depth of the eclipse curve and 
vary its central phase, but leave its shape (duration, limb- 
crossing time) unchanged. At all wavelengths, we fit a 
baseline curve plus the eclipse curve simultaneously via 
multi- variable linear regression, but the nature of the base- 
line curve varies with wavelength (see below). Since the 
linear regressions cannot fit a variable central phase, we 
perform the regressions separately for each of a series of 



central phase values, and we pick the solution having the 
minimum % 2 . 

We found that a linear baseline was adequate for 
the eclipse fits at all wavelen gths except 8.0 /im where 
the well-known detector ramp ( Charbo nneau et all 120051: 
Deming et al.ll2006t Harrington et al.ll2007t Knutson et alj 



2008. 120091: iDesert et al.ll2009h exhibits a quasi-logarithmic 
shape. We model the 8 /im ramp, i?(^>), as a sum of a linear 
and logarithmic term in phase (cf>): 



R{4>) = a + a\4> + a 2 ln(4> - <j) ) 



(1) 



We adopt multiple values for </>0j and solve for ao, a\, 
and ct2 by linear regression at each adopted value for (j>Q. 
Note that 0o is a phase offset used to facilitate the ramp fit, 
and is not related to the orbit of the planet. We choose 
the best fit from the 2-D grid of 4>q and eclipse central 
phase values based on the minimum % 2 . We found that 
this model of the ramp provides consistently excellent fits, 
but the best-fit ai and 0o can be degenerate in the sense 
that different combinations can produce indistinguishable 
ramps. Fortunately, we did not find the eclipse depth and 
central phase to exhibit significant degeneracies with the 
ramp parameters. Figure 4 shows the 5.8 and 8/im data 
before ramp removals, with the best-fit ramps and eclipses 
overlaid. Figure 5 shows the eclipse fits at 5.8 and 8/im, 
in comparison to binned data with the baseline and ramp 
effects removed. 

Note that a possible temporal drift in intensity at 3.6 
and 4.5 /im is a phenomenon physically distinct from the 
pixel position effect. After correcting intensity for pixel po- 
sition, we include a linear baseline when fitting the eclipse 
curves. We have added these linear baselines to the total 
decorrelation function illustrated in Figure 2. The slope 
of the baseline at 3. 6 /mi is 0.018% ± 0.003% per hour. 
iKnutson et alj (|2009f ) found a linear increase at 3.6 /im of 
similar magnitude. We found that all of the intensity vari- 
ations at 4.5 /im were fully accounted for by changes in 
pixel position, and the temporal term was not significant. 

Best-fit eclipse depths and errors, in units of the stellar 
flux (contrast units) are given in Table 1. 

3.4. Error Estimation 

We esti mate the errors usin g the bootstrap Monte Carlo 
technique (|Press et al.l 1 1992T ) . The bootstrap technique 
generates synthetic data sets using the residuals from the 
best-fit model, and permutes them to make new data. 
Each new boostrap data set is constructed as follows. To 
construct N new data points, we start with the i — 1, N 
points from the best-fit curve. We draw a residual ran- 
domly from the pool of original residuals, add it to the 
i th best-fit curve point, return that residual to the pool, 
and draw again until we have created a data set of N new 
points. We make 50,000 data sets using this procedure. 

For each new bootstrap data set, we repeated the entire 
fitting procedure, with the exception of the pixel position 
fits. We did this for the 50,000 bootstrap data sets at 
each wavelength, and tabulated histograms of the eclipse 
depth, baseline parameters, and central phase. These his- 
tograms are very close to Gaussian distributions, and their 
standard deviations give estimates of the error in the fit 
parameters. 
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We also estimated errors using the res idual permutation 
method described by ISouthwortbl (|2008l ). This method - 
sometimes called the 'prayer bead' method - is similar to 
the bootstrap technique, except that it preserves the or- 
der of the residuals, and is therefore more sensitive to the 
presence of red noise. In most cases the error from the 
permutation method was very close to the bootstrap er- 
ror. The adopted error for a given parameter was taken to 
be the greater of the values from the bootstrap and resid- 
ual permutation method, and Table 1 includes these errors 
for the eclipse depth and central phase. 

4. RESULTS AND DISCUSSION 
4.1. Atmospheric Temperature Structure 

Figure 6 shows our contrast valu es plotted versu s wave- 
length, and also the results from iKnutson et al.l {2008) 
for HD 209458b, the archetype of an inverted atmosphere. 
We include the contrast predicted b y two 1-D models of 
the H AT-P-lb planetary atmosphere (|Fortnev et al.|[2005l 
2008]), both invoking re-emission of stellar irradiance on 
the dayside only. The solid line is the nominal model, with- 
out a temperature inversion. This nominal model is self- 
consistent and non-gray, with solar metallicity, and was 
calculated as in lFortnev et al. ( 2008). It has a bond albedo 
of 0.067 and a dayside effective temperature of 1512K. This 
model has no inversion because we find that it is too cool 
to allow gas-phase TiO. 

The dashed line is a weakly inverted model, also of so- 
lar metallicity. It is not self-consistent; it has a constant, 
ad-hoc temperature inversion of dT/d log P = — 30K. How- 
ever, this model is constrained to have the same effective 
temperature (1512K) as the non-inverted model. Since the 
temperature gradient of the inverted model is shallow, the 
emission features are weak, and it falls close to the contrast 
expected from a 1500K blackbody (dot-dashed line). In 
all instances , we u sed the planetary and stellar radii from 
IWinn et all {20071) . and a Kur iicz 6000/4.5/0.0 m odel to 
represent the stellar spectrum ([Torres et al.ll200'8l) . 

The rotation of HAT-P-lb should be tidally locked to 
its orbit, even though its orbital period is longer than 
fo r many planets in t he hot Jupiter class. Using Eq.(l) 
of lGuillot etaLI ([1996), we calculate a spin-down time of 
8 x 10 6 years, starting with Jupiter's rotation rate and 
adopting a very conservative Q-value (10 6 ). Since this 
spin-down time is much less than the age of the system 
([Torres et al.ll2008() . we expect tidal locking of the planet's 
rotation. 

Based on tidally-locked rotation, the maximum dayside 
temperature of HAT-P-lb is 1550K, assuming zero albedo, 
a uniform temperature over the dayside hemisphere, and 
no transport to the nightside. This is too cool to produce 
an inversion using TiO / VO absorption, so we first investi- 
gate the non-inverted model. Emission at the three longest 
IRAC wavelengths can arise from levels hi gher in the at- 
mosp here than does the 3.6 /im radiation (|Burrows et alj 
2007J), and these three channels exhibit a contrast that is 
moderately higher than the non-in verted model. Com- 
paring to the HD 209458b results (jKnutson et alj l2008f) 
shows good agreement at 3.6 and 8.0 /im, but HAT-P-lb 
exhibits a contrast at 4.5 and 5.8 /im that is intermediate 
between HD 209458b and the non-inverted model. This 
seems qualitatively consistent with a moderate inversion, 



perhaps produced by an absorber other than TiO/VO. 

We integrated the flux from each planetary atmospheric 
model and the stellar model atmosphere over the IRAC 
bandpass functions. Dividing these integrated fluxes at 
each wavelength produces expectation values for the ob- 
servations. (These are not illustrated on Figure 6, but 
they fall very close to the contrast values at the band- 
center wavelengths.) We used these expectation values to 
calculate the x 2 value for the observations compared to 
each model. For the inverted model (dashed line), \ 2 1S 
3.9, whereas it is 11.0 for the non-inverted model. A value 
as high as 11.0 will occur only 2.6% of the time if the non- 
inverted model is an accurate description of HAT-P-lb's 
atmosphere. This level of confidence is not sufficient to 
rigorously eliminate the non-inverted model, but it does 
indicate that the inverted model is preferable. 

Within the errors, our results for HAT-P-lb can also 
be reproduced using a blackbody spectrum for the planet 
(dot-dashed line on Figure 6, x 2 = 3.5). 

The luminosity of HAT-P-lA is 1.48 times solar 
(|Torres et al.l l2008h . and it receives a stellar flux ~ 2/3 
of the HD 209458b case. If HAT-P-lb absorbs with zero 
albedo and re-emits uniformly but only on the dayside 
hemisphere, then we expect a maximum temperature of 
1550K. If the planet's emission is uniform over both hemi- 
spheres, we expect an observed temperature of 1300K. Our 
observations at secondary eclipse are best described by a 
blackbody having a temperature of 1500 ± 100K, where 
we have factored in the random error of our observations 
as well as uncertainty in the stellar and planetary pa- 
rameters. However, m ost of the flux is probably emitte d 
at shorter wavelengths (|Barman 2008; Seager et al.ll2005l ). 
not directly observable using Spitzer. At face value our re- 
sults suggest that redistribution of the stellar irradiation 
by dynamics may be inefficient for this planet. 

4.2. Orbital Eccentricity 

HAT-P-lb orbits a star in a wide visual binary, and 
this circumstance can have significant consequences for 
the orbital dynamics of the planet. The inclination and 
eccentricity (e) of the planet's orbit can in principle un- 
dergo oscilla tions and long term evolution due to the Kozai 
mechanism {Fabrvckv fc Tremainel 120071 ). Misalignments 
between the planet's orbital inclination and the rotation 
axis of the star can result, a nd can be observed using 
the Rossiter-McLaughlin effect. iJohnson et alj {2008) find 
that the angle between the sky projections of stellar spin 
axis and the orbit normal for HAT-P-lb is 3.7 ± 2.1 de- 
grees. They also put an upper limit on the eccentricity of 
0.067 with 99% confidence. These values suggest that any 
Kozai oscillations of HAT-P-lb have largely damped out. 
Our results for the secondary eclipse timing will further 
strengthen that conclusion, as discussed below. 

The timing of the secondary eclipse is exquisitely sen- 
sitive to the eccentricity of the orbit, with one ambiguity 
being the value of u>, the argument of periastron. When 
our line of sight aligns with the minor axis of the planet's 
orbit (u> = or n), then the secondary eclipse will not be 
centered on phase 0.5 unless the orbit is circular. When 
our line of sight aligns with the major axis of the orbit 
(u> = 7r/2 or 3tt/2), departures from circularity will affect 
the duration of the eclipse, but not the central phase. The 
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phase of s econdary eclipse therefor e constrains the value 
of ecosw (|Charbonneau et aT1l2005h . 

We observed two eclipses, each at two wavelengths si- 
multaneously. The difference in best fit phases for the 
same eclipse observed at different wavelengths is approxi- 
mately consistent with our errors (Table 1). The larger dif- 
ference occurs at 3.6 and 5.8 /mi, where the central phase 
difference is 0.5016 - 0.4992 = 0.0027. Because the noise 
at each wavelength is independent, the error on the differ- 
ence in the two phases is the quadrature sum (= 0.0014) of 
the phase errors at the two wavelengths. Hence the phase 
difference is less than a 2a variation. 

Weighting the central phase of the eclipse at each wave- 
length (Table 1) by the inverse of its variance, we find a 
mean value of 0.4999 ±0.0005. Considering the 55 seconds 



of light travel time across the orbit, we expect to find the 
eclipse at phase 0.500014. Accounting for the light travel 
time and the error in the observed eclipse phase, we derive 
a 3cr upper limit of |e cos u>\ < 0.002. If this planet's or- 
bital eccentricity was affected by Kozai oscillations in the 
past, they have damped to the point where the orbit is 
closely circular. 

This work is based on observations made with the 
Spitzer Space Telescope, which is operated by the Jet 
Propulsion Laboratory, California Institute of Technology, 
under a contract with NASA. Support for this work was 
provided by NASA. We are grateful to the anonymous ref- 
eree for thoughtful comments that improved this paper. 

Facilities: Spitzer. 
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Fig. 1. — Upper panel: Photometric intensity versus Y-pixel position for 3.6 /xm photometry, before correction. Pixel centers are at integer 
coordinate values. The red points represent lower intensities due to approaching saturation for this bright star as it moves closer to pixel 
center in the Y coordinate, and those points were not used in our analysis. The blue line is the fitted correction function. The X-position of 
the star (not illustrated) was approximately pixel 79. (Both X- and Y-pixel coordinates are 1-based.) Lower panel: Intensity versus orbital 
phase for the 3.6 fim photometry, before the correction. The red points are not used in the analysis. Note that the eclipse is visible in this 
plot, but see Figs. 2 and 3 for a clearer view. 
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Phase 

Fig. 2. — Photometry at 3.6 fim (upper panel) and 4.5 fim (lower panel), before decorrelation with image position. For clarity, each plotted 
point represents the binned average of the original photometry. In the upper panel each bin represents 10 exposures (about 2.4 minutes per 
bin). In the lower panel, each bin represents 20 exposures (about 2.2 minutes per bin). The lines show the decorrelation functions. Note that 
these functions are smooth when plotted as a function of pixel phase, but because of pointing jitter they exhibit fluctuations when plotted 
here as a function of orbital phase. The eclipses are clearly visible near orbital phase 0.5, as differences between the points and lines. The 
decorrelation function at 3.6 /im includes a linear drift as a function of time (see text), but all of the variations at 4.5 fim can be attributed 
to changes in pixel position. 
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Fig. 3. — Binned photometry at 3.6/jm (upper panel, 10 exposures per bin) and 4.5 fim (lower panel, 20 exposures per bin), corrected for 
the pixel phase effect (see Figures 1 &; 2, and text). The solid red line shows the fitted secondary eclipse, and the dashed red lines show the 
± 1<t range on the eclipse depth, from the bootstrap Monte-Carlo trials (see Table 1, and text). The error bars give the standard deviation 
of the mean, based on the internal scatter within each bin. 



Spitzer Photometry of HAT-P-lb 



9 



0.48 0.49 0.50 0.51 0.52 
••* y: 




. . i •■ ■ ■ i i i i . . 

0.48 0.49 0.50 0.51 0.52 

Phase 



Fig. 4. — Photometry at 5.8 (upper panel) and 8.0 (im (lower panel) before removal of the detector ramp. The points are unbinned. The solid 
red lines are the best fit baseline ramps plus eclipse curves, obtained via linear regression (see text). A linear ramp was used at 5.8 /im, whereas 
at 8.0 fj,m the ramp is comprised of a linear plus logarithmic term, and the best-fit (see text) is 1.005 X 10 — 2 ln(if> — 0.465) — 0.07</> + 1.07063, 
where tf> is phase. 
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Fig. 5. — Binned photometry at 5.8 /im (upper panel) and 8 /an (lower panel), with the baselines removed. The bin size is 16 points at 
both wavelengths, about 3.6 minutes per bin, but the eclipse fit used the unbinned data (sec Figure 4). The solid red lines show the fitted 
secondary eclipse, and the dashed red lines show the ± lc range on the eclipse depth, from the bootstrap Monte-Carlo trials (see Table 1, 
and text). The error bars give the standard deviation of the mean, based on the internal scatter within each bin. 
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Fig. 6. — Contrast values for HAT-P-lb (solid circles) versus wavelength. The open squares are for HD 209458b, from Knu tson et all l(200Sly 
The two obser ved points at 3.6 /im have been offset slightly in wavelength to prevent overlap of their error bars. The solid line is a model 
for the planet IIFortnev et al. 2008) with no temperature inversion, and with re-distribution of stellar irradiance over the dayside only - with 
none to the nightside. The dashed line is for a model with a modest temperature inversion (see text). In the inverted model, very modest 
decreases in contrast can be seen at the wavelengths of contrast peaks in the non-inverted m odel. The dot-dashed line is the contrast for a 
1500K blackbody having the same radius as HAT-P-lb I Winn et al.ll2007l ; iTorres et al.ll2008T) . 
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Table 1 

Fitted Eclipse Depth, Central Phase and Heliocentric Julian Date (HJD) of Mid-eclipse. 



Wavelength Eclipse Depth Central Phase HJD-2454000 

3.6 0.080% ± 0.008% 0.5016 ± 0.0008 464.4228 ± 0.0036 

4.5 0.135% ± 0.022% 0.4991 ± 0.0010 102.7229 ± 0.0045 

5.8 0.203% ±0.031% 0.4992 ± 0.0012 464.4121 ± 0.0054 

8.0 Lim 0.238% ±0.040% 0.4986 ± 0.0010 102.7206 ± 0.0038 



